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We consider a single vortex in a superfluid Fermi gas in the BCS-BEC crossover regime near a 
Feshbach resonance. The effect of the molecular Bose-Einstein condensate upon the vortex struc- 
ture is discussed within the mean field approximation at zero temperature. Using the self-consistent 
Bogoliubov-de Gennes equation of the fermion-boson coupled model, we calculate density distribu- 
tions of atoms and molecules. As the number of the molecules increases, both atomic and molecular 
density changes from BCS-like distribution to BEC-like. We also study the change of the vortex 
core size in the crossover regime. 
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I. INTRODUCTION 

The evidence of fermionic superfluidity in trapped 
atomic gases has been clearly shown in the recent ex- 
periments P, 0, 0, S • The significant feature of these 
systems is that the strength of the inter-atomic interac- 
tion is manipulated via a Feshbach resonance (FR) phe- 
nomena. The crossover, which has been discussed sev- 
eral decades Q, |^ , between a weak coupling Bardeen- 
Cooper-Schrieffer (BCS) superfluid and a Bose-Einstein 
condensate (BEC) of pre-formed pairs can be realized in 
these systems. 

The FR phenomena involve the scattering of atoms 
from open channel states into a molecular bound state 
formed from neighboring closed channel states. The so- 
called "pre-formed pair" in the conventional BCS-BEC 
crossover theory is here equivalent to the substantial 
molecule in the closed channel, which has been observed 
as a BEC on the one side of the FR. 

In this paper we discuss the property of a single vortex 
in the BCS-BEC crossover regime at zero temperature. 
Since the vortex states in BCS superfluids and in BECs 
are significantly differ in their distributions of particles, 
observing a vortex state in a Fermi gas provides the in- 
formation about superfluidity in the crossover regime. In 
the case of a vortex in the BCS state, the energy gap is 
suppressed at the vortex core, and therefore the particle- 
hole symmetrical modes, which are localized at the vor- 
tex core, exist in the vicinity of the Fermi surface 0. 
Therefore, the atomic density is finite at the center of 
the vortex core. On the other hand, the particle density 
in a BEC vanishes at the vortex core, directly reflecting 
a singularity of the order parameter. 

We are also interested in the vortex core size in the 
crossover regime. The core size in a BCS superfluid is 
given by ^bcs ~ hvp/ A, where vp is the Fermi velocity 
and A is the energy gap. As the attractive interaction 
becomes strong, the growth of the energy gap makes the 
core size smaller for a molecular BEC, the core 

size is written by ^bec ~ VvS^rrtMOM, where hm is the 



number density of molecules and aM is the s-wave scat- 
tering length of a molecule, which diverges at resonance 
as well as that of an atom 17, 15] . So ^bec a-lso becomes 
small near the resonance. Then, how does the core size 
change in the crossover regime? 

There are several theoretical papers concerning vortex 
states in Fermi gases JJ, .14, ,15, 16, 17J. These 
papers, however, deal with only Fermi gases of which 
interaction strength is manipulated, and the effect of a 
molecular BEC is not considered. In this paper, we start 
with an atom- molecule coupled model 0, figl l20l | , and 
calculate density distributions and order parameters of 
both atoms and molecules. We use the mean field ap- 
proximation to simply deal with the effect of a molecular 
BEC. 



II. FORMALISM 

We consider atomic gases of two atomic hyperfine 
states (labeled by cr =t|), which are coupled to a molec- 
ular two-particle bound state. The energy of a bare 
molecule relative to that of two bare atoms is denoted by 
2j/. The model Hamiltonian of the fermion-boson cou- 
pled system is given by 



H 



(7 

+ 0^(r) {7^M(r) +2i/}0(r) 



(1) 



where tpa and <f> are the field operators of atoms and 
molecules, respectively. The first and second terms in the 
integration represent the bare Hamiltonian of atoms and 
molecules, respectively, where 7iA(r) = — (ft^/2m)V^ -f- 
Va{y), HM(r) = -{h? /Am)\/'^ + Vm{y), m is the atomic 
mass and Va,m{y) are the trapping potential for atoms 
and molecules, respectively. The last term in Eq. Q rep- 
resents the atom-molecule coupling associated with FR. 
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Its contribution to the effective atom-atom interaction is 
—g'^/2v. This means that the interaction is manipulated 
by changing u. So we investigate the u dependence of 
this system, especially at around v ^ {). Although there 
exist non-resonant inter-atomic and intermolecular inter- 
actions, we assume that the effect of the resonant process 
is so dominant that the other processes are negligible. 

To conserve the total number of particles, we use a 
single chemical potential /i and work with H = H — fiN, 
where N is the number operator given by 



N 



dr 



(2) 



Here, we introduce the mean-field order parameters 
■p(r) = (■0|(r)'0|(r)) and 0(r) = {(j){r)), which cor- 
respond to the Cooper pair amplitude and molecu- 
lar order parameter, respectively. The number den- 
sity of molecules are given by riA/(r) = |(/)(r)p, as 
usual. Though we neglect the bare atom-atom interac- 
tion, the effective interaction via FR supports the pairing 
of atoms. Since <j>{r) is related to an energy gap as we 
will show below, the pair amplitude remains finite as long 
as a molecular BEC exists. The relation between these 
order parameters is given by the equilibrium condition 
ih{d(t){r) / dt) = ([0(r),7i]) = 0, from which it follows 
that 



[Ha/ (r) + 2i^- 2^] cf>{r) + gV{v) = 0. 



(3) 



As in the case of the well-known mean field BCS theory, 
the "off-diagonal energy" should be defined as the energy 
gap, i.e., Eq. implies that Aoff (r) = g(t){v). Then, the 
conventional Bogoliubov-de Gennes (BdG) approach is 
used in a similar manner: the field operator of the atom 
is transformed with 



V'](r) 



E 



(4) 



and the BdG equation is obtained by 









■ u.j{r) ' 



















.(5) 



In this notation, the pair amplitude and the density at 
zero temperature are written by V{v) = — it*(r)wj(r) 
and n^(r) — |wj(r)p, respectively. To avoid the ul- 
traviolet divergence in 7^(r), we introduce a cutoff energy 
LJc which is in the order of the Fermi energy, and redefine 

Then, all we have to do is to solve Eqs. lO and (O 
self-consistently. The chemical potential is determined 
so that the total number of atoms A'tot = Na + 2Nm is 
conserved, where Na,m = J c^rn^.Af (r). 



III. VORTEX STATE IN THE CROSSOVER 
REGIME 

We assume a cylindrical optical trap and approximate 
it with a two-dimensional harmonic potential: Vyi(r) = 
(l/2)m[j2(a;2 + y^) and VA,f(r) = 2VA(r). For simplic- 
ity, we neglect the z dependence of the order param- 
eters. When considering the microscopic vortex struc- 
ture, however, the degree of freedom in z direction is 
important. So we assume that gases are confined in the 
length Lz, and impose the periodic boundary condition. 
We have numerically confirmed that the results does not 
depend on but on Ntot/Lz when Lz > oho, where 
flfjo = \/ h/2mu is the harmonic oscillator length. 

A single particle state of a bare atom in this potential 
is given by e^'=^e*'^(r/aHo)'>ci!' (rV2a|o)e-'''/4"Ho hav- 
ing an energy eigenvalue ekin = h^{2n + l + k^a^fj), 
where (r, 6*, z) is the cylindrical coordinate and /^n{x) is 
the generalized Laguerre polynomial function. The in- 
dices specifying the energy are given by ^ = 0, ±1, ±2, • • • , 
n = 0,1,2,---, and fcano = (27raHo/-^z)'^fc = 
where Uk = 0,±1,±2, • • • . By counting the number of 
the eigenstates below ep, the Fermi energy as a func- 
tion of the number of atoms is given by e_F(AOi) = 
huj{15koNA/l6)'^^^ ■ The characteristic energy of the sys- 
tem is defined as Ep = epiNtot), which corresponds to 
the Fermi energy in the BCS limit. 

For a vortex state in equilibrium, we assume the form 



Wfe;n(r) = e" 



ikz ilO /\ 



(6) 
(7) 
(8) 



where Ukin{j'),Vkin{f) and (/)(r) are real functions, and 
n = 0, 1,2,--- is the radial quantum number. In the 
BEC limit, by substituting gV = Q in Eq. © the lowest 
energy state of a molecular BEC with a single vortex 
is given by (j){v) = Nm I'^'^a^QLze''^ fo{r / a-ao), where 
/o(r) = ^/2re^^ 1"^ . Here, the chemical potential is also 
determined by Eq. Q as n = v + Huj. 

We consider a gas of ^Li atoms in a trap with uj — 
27r X 300 Hz and ajjo = 1-7 /im. The line density are 
set as Ntot/Lz = 30 /zm~^ and 200 /Ltm~-^, leading to 
Ep — 9.7hijj and 21Hu}, respectively. Though the particle 
density in our calculation is about 10 times smaller than 
that in experiments, Ep/hw are in the same order. As 
for the coupling constant, we consider a narrow FR and 
set g\/n — Ep and ZEp^ where n is the mean density of 
total atoms: n — Ntot/2Tra'^QLz- 

Figure ^a) shows the ly dependence of the chemical 
potential. Since we neglect the bare atom-atom and 
molecule-molecule interactions, the trapped gas behaves 
as an ideal gas far from the resonance in each side as 
shown in Fig. Qfa), i.e., in the BCS hmit ^ goes to Ep, 
while ^ ^ V + hui in the BEC limit. In Fig. ^b), the 
numbers of atoms and molecules are plotted as a function 
of v. It is clearly shown that where the crossover occurs 
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FIG. 1: (Color online) (a) The v dependence of the chem- 
ical potential for Ef ~ 9.7huj. The inset shows the results 
for Ef ~ 21huj. The chemical potential far from the reso- 
nance goes to the asymptotic forms; /i = Ef in the BCS limit 
{ly > 0) and = ly + huj in the BEC limit {u < 0). (b) The 
V dependences of A'^^ (thick lines) and 2Nm (thin lines) . The 
total density is always conserved: Na + 2Nm = Ntot- The 
intersections of Na and 2Nm depend upon g^/n/Ep; two in- 
tersections in the case of gy/n/EF = 3 are coincide with each 
other. 

in the parameter space v/Ep is determined by the inter- 
action energy g^/ri relative to Ep, while the width of the 
crossover region depends on the couphng constant g. 

IV. CORE STRUCTURE 

Figure 13 shows the density distributions of (a) atoms 
and (b) molecules in the case of Ep — 21huj and .gy^ = 
3Ep. Both profiles are normalized by the correspond- 
ing numbers of each particles. We also plot the density 
distribution of a non- interacting BEC, |/o('')P, in both 
figures. As we mentioned above, the atomic distribution 
far from resonance in BCS side {v/Ep = 1.5) has no 
explicit hole associated with a vortex. The broad distri- 
bution in radial direction is also the feature of a fermionic 
profile. As v decreases, however, these fermionic features 
disappear and the atomic distribution becomes BEC-like, 
i.e., the cloud becomes narrower in radial direction and, 
moreover, the hole at the vortex core clearly appears. 
This result is consistent with that of the one-channel 



FIG. 2: (Color online) The density distributions of (a) atoms 
and (b) molecules in the case of Ef = 21hu) and g^/n = 3i?F. 
The density profile of an ideal Rose gas is also plotted in 
both figures. Both atoms and molecules change from broad 
distribution at v/Ef = 1.5 to narrow one at v/Ef ~ —1.0. 
The fraction of atoms localized at the core decreases as u de- 
creases, while the molecular density at a vortex core is always 
zero. 



model [13, which argues that the density decreases at 
the vortex core when a BCS superfluid goes into the BEC 
regime. Our calculation shows that this density depletion 
does not merely come from the increase of the molecu- 
lar fraction. The atomic profile itself also turns into a 
BEC-like distribution. 

The existence of a molecular BEC leads to these 
changes in atomic distribution. Since the molecular den- 
sity always vanishes at vortex core as shown in Fig.[21[b), 
the atomic density is also suppressed there when the 
number of molecules becomes large. On the other hand, 
the molecular distributions in the BCS side are also 
strongly affected by atoms, being broader than those in 
BEC side. In the crossover regime, atoms and molecules 
are strongly coupled to each other, and both gradually 
turn from BCS-like to BEC-like distribution. 

Next, we discuss the change of order parameters. Fig- 
ure |31 shows the V dependence of (a) the pair amplitude 
7^(r) and (b) the molecular order parameter 0(r), or the 
effective energy gap. These are also the results of the 
calculation with Ep — 21/iw and g^/n = 3Ep. As v 
decreases, the molecular order parameter monotonically 
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FIG. 3: (Color online) The spatial distribution of the mean 
field order parameters, (a) The pair amplitude of atoms 
changes the v dependence at v jEp ^ 1.0, while (b) the molec- 
ular order parameter monotonically grows as v increases. 

grows and becomes narrower as well as the number den- 
sity. On the other hand, the fermionic pair amplitude 
becomes largest at vjEp = 1.0 and goes to zero in both 
BCS and BEG limit. As v decreases in the BCS side, the 
pair amplitude increases since the effective energy gap, or 
(^(r), increases. In the BEG side, however, an atomic pair 
turns into a molecule and the pair amplitude becomes 
smaller as the number of molecules increases. The v de- 
pendence of the pair amplitude changes at v jEp = 1.0, 
where the energy gap is equal to the chemical potential. 

In Figs.|2Ia) and (b), the order parameters have max- 
ima at a nearly same radius. This means that the core 
size does not strongly depend on v. In the BGS side, 
the core region of V(r) becomes smaller as v decreases 
(see, vjEp = 1.5 ~ 1.0 in Fig. (SJa)). This fact is con- 



sistent with the BGS theory. When the line density of 
total atoms is set larger as in the experiments, leading 
to a broader distribution in radial direction, this change 
must be clear. In the region where v/Ep < 1.0, how- 
ever, the core size does not become smaller any more. 
The core size in this region is determined by the scale of 
spatial variations of the single-particle wavefunction with 
the energy Ep, which is in the order of ajjo- In the case 
of free atoms, this size corresponds to the Fermi wave 
length 1/kp = h/y/2mEp. 

As for the molecular order parameter, the core region 
of (j){r) in BEG limit corresponds to the harmonic oscil- 
lator length since we neglect the bare molecule-molecule 
interaction. In the crossover regime, the atomic pair is 
expected to affect the core size of the molecular order pa- 
rameter. It seems, however, that the effect of the change 
in the V{r) on the core size in (/)(r) is very little in our 
results. 



V. CONCLUSION 

We have studied a vortex state in a superfluid Fermi 
gas near a Feshbach resonance, using the atom-molecule 
coupled model and the mean field approximation. Gou- 
pled to each other, both atomic and molecular density 
distributions change from BGS-like to BEG-like in the 
crossover regime. The depletion of the total particle den- 
sity in BEG side does not merely come from the increase 
of the number of molecules. Affected by a molecular 
BEG, a hole at the core in atomic density clearly ap- 
pears. We have also studied the change of a vortex core 
size. In the BGS side, the change of the core size is con- 
sistent with the BGS theory. But becouse of the trapping 
potential, the core size dose not show drastic change in 
the crossover regime. 
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